Collisional energy loss of a fast heavy quark in a quark-gluon plasma 
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We discuss the average collisional energy loss dE/dx of a heavy quark crossing a quark-gluon 
plasma, in the limit of high quark energy E 3> M 2 /T, where M is the quark mass and T <C M 
is the plasma temperature. In the fixed coupling approximation, at leading order dE/dx oc a 2 , 
with a coefficient which is logarithmically enhanced. The soft logarithm arising from t-channel 
scattering off thermal partons is well-known, but a collinear logarithm from u-channel exchange had 
previously been overlooked. We also determine the constant beyond those leading logarithms. We 
then generalize our calculation of dE/dx to the case of running coupling. We estimate the remaining 
theoretical uncertainty of dE/dx, which turns out to be quite large under RHIC conditions. Finally, 
we point out an approximate relation between dE/dx and the QCD Debye mass, from which we 
derive an upper bound to dE/dx for all quark energies. 

PACS numbers: 12.38.Bx 



I. INTRODUCTION 

Jet quenching, as first anticipated by Bjorken [1], is 
a crucial probe of the state of matter created at the 
Relativistic Heavy Ion Collider (RHIC). The quenching 
of light hadron spectra at large p± [2, 3] is usually at- 
tributed to the radiative energy loss of the (light) parent 
parton when crossing the hot or dense medium. On the 
other hand, recent data on heavy flavour quenching [4, 5] 
suggest that for heavy quarks, purely radiative energy 
loss might be insufficient to explain the observed atten- 
uation. This has renewed the interest in the collisional 
part AE co ii of the parton energy loss, in particular in 
the case of a heavy quark [6] . It is not clear yet whether 
such a collisional contribution can help understanding the 
data on heavy flavour quenching. For instance, although 
some studies support the possibility of a quite large colli- 
sional loss [7-9] , a recent comparison between collisional 
and radiative losses [10] indicates that the collisional con- 
tribution to parton energy loss might be small (~ 20%) 
compared to the radiative one, even for heavy quarks. It 
should however be noted that a relatively small average 
collisional loss might be compatible with an important ef- 
fect of collisions on quenching [11, 12], due to differently 
behaved collisional and radiative energy loss probability 
distributions. 

In any case, in the present (unclear) situation where 
the importance of collisional energy loss is reconsidered, 
we believe it is not useless to state the correct result for 
the (average) collisional loss of an energetic heavy quark. 

A basic quantity required to study collisional quench- 
ing is the rate of energy loss per unit distance, dE/dx, 
of a parton produced in the remote past and travelling 
in a large size medium, as studied in Refs. [1, 13, 14]. 
For heavy ion collisions, where a parton initially pro- 
duced in a hard subprocess crosses a medium of fi- 
nite size L, we expect deviations from the linear law 
AE co u(L) = (dE/dx) ■ L [15]. Nonetheless, knowing the 



'asymptotic' rate dE/dx is a prerequisite before attempt- 
ing a more refined evaluation of AE co u , including in par- 
ticular finite size effects. 

So far, the most detailed calculation of dE/dx for a 
heavy quark in a quark-gluon plasma (QGP) is due to 
Braaten and Thoma (BT) [14], and is based on their 
previous QED calculation of dE/dx for a muon [16]. As 
we recently analyzed [17], the BT QED calculation re- 
lies on an assumption on the momentum exchange q in 
elastic scattering, namely q <C E, which is incorrect in 
the domain E > M 2 /T. Therefore, the BT results for 
dE/dx, both in QED [16] and in QCD [14], need to be 
corrected in this limit. This was done in [17] in the QED 
case. Here we consider the QCD case of a fast heavy 
quark. 

We assume the heavy quark to be produced in the hard 
partonic subprocess of the heavy-ion collision, i. e. to be 
present in the 'initial' stage of the QGP evolution. It 
then crosses the QGP (supposed thermally equilibrated) 
by losing some energy, before hadronizing (within a jet) 
into a heavy-flavoured hadron. The situation we consider 
should be appropriate when dealing with heavy-quark 
tagged jets, where the energy loss is defined as the differ- 
ence between the initial and final energies of the flagged 
heavy quark. This is different from light parton (or un- 
tagged heavy quark) energy loss, which should be defined 
(at the partonic level) as the energy difference between 
the leading partons in the initial and final states. Indeed, 
without tagging, the flavour or even the type (quark or 
gluon) of the initial energetic parton does not have to be 
conserved in the 'energy loss' process. 

In section II we study the fixed coupling approximation 
which allows us to closely follow the lines of our QED cal- 
culation [17]. In section III we show how the result with 
a running coupling can be simply inferred, and discuss 
the theoretical uncertainty. The latter is illustrated by a 
numerical estimate in section IV. 
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FIG. 1: Amplitudes for heavy muon (quark) elastic scattering in a QED (QCD) plasma. A curly line denotes a photon (QED) 
or a gluon (QCD). The amplitude (d) is specific to the QCD case. The blob in (a) and (d) denotes the resummed hard thermal 
loop boson propagator, which is necessary to screen the t-channel contribution in the infrared. 



II. FIXED COUPLING APPROXIMATION 

In this section we assume a fixed coupling, which is jus- 
tified when the running coupling does not change much in 
the range of probed momenta (we will specify this con- 
dition below). We calculate, in the limit E > M 2 /T, 
the heavy quark collisional loss dE/dx beyond logarith- 
mic accuracy. For the sake of clarity we first focus on the 
leading logarithms before calculating the constant next 
to those logarithms. 

A. Leading logarithmic terms 

The leading logarithms in dE/dx for a heavy quark in 
the quark-gluon plasma have the same origin as in the 



QED analog of a muon crossing an e 7 plasma. We thus 
make a detour to QED and outline the derivation of [17], 
which is straightforward at leading logarithmic accuracy. 
The QCD result for dE/dx of a heavy quark will then be 
obtained from simple considerations. 



1. Muon energy loss in a QED plasma 

At leading order the energy loss of a muon (of mass M 
and energy E) arises from elastic scattering off thermal 
electrons or positrons (Fig. la) and photons (Figs, lb- 
c) . Scattering off electrons corresponds to i-channel ex- 
change, whereas scattering off photons (Compton scat- 
tering) receives contributions from s and u-channels. The 
muon energy loss is given by [16] 



dE^ 
dx 
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(2.1) 
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where oj = E — E' ' . The tree- level amplitude Mi corresponds to scattering off a thermal particle of type i 
Each |.Mi| 2 is summed over initial and final spin states, and we divide by the degeneracy factor d = 2 of the incoming 
muon. Furthermore, nj(fc) = (exp(fc/T) =p 1) _1 is the thermal distribution of the target particles, and fn = 1 ± rij 
accounts for the Bose enhancement or Pauli blocking for the scattered state. We also use the shorthand notation 
f k = fd^/(2nr. | 



In the E > T limit, (2.1) can be simplified to [17] 



dx 



n,(fc) 
2k 



dt (-t) 



doj_ 
~dt 



(2.2) 



where di is the degeneracy factor of the target particles, 
and 



dai 
~dt 



1 
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is the corresponding differential cross section. We also 
define s = s — M 2 , as well as for later reference u = u — 
M 2 . In Eq. (2.2) the bounds on t are set by kinematics, 
imin = —s 2 /s, and i max = 0. We will focus on the limit 
E > M 2 /T, which implies s = (P + K) 2 = M 2 + 2PK ~ 
O (ET) > M 2 , so that i min ~ -s. 



To obtain the leading logarithm from the t-channcl 
contribution we can assume \t\ <§C s and approximate 

]T ^|M i | 2 ^32e 4 ^. (2.4) 

i=e± spins 



This contributes to (2.2) as 
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To obtain the r.h.s., we replaced i m j n ~ — s — > —ET, 
which is justified to logarithmic accuracy. Strictly speak- 
ing, with i max = the integral in (2.5) would be infrared 
divergent. As is well-known, this divergence is screened 
by medium effects. To logarithmic accuracy it is suffi- 
cient to take i max = — m D as an effective infrared cut-off, 
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where bid = el is the Debye screening mass in the 
QED plasma [18]. The logarithm in (2.5) thus arises from 
the domain m 2 D <C \t\ <C s ~ ET. A more accurate de- 
scription of screening requires using the resummed hard 
thermal loop (HTL) [18, 19] photon propagator in the 
t-channel (instead of the bare one with effective cut-off), 
as pictured in Fig. la. This is needed to control the con- 
stant ~ O (1) in (2.5), as will be recalled in section II B. 
Finally, in order to obtain (2.5) we used 



2k 



48 



(2.6) 



The contribution from Compton scattering to the 
muon energy loss brings another logarithm, arising from 
the square of the w-channel amplitude, more specifically 

[17]. In this do- 



from the domain u r 



< u < u r 



main the Compton scattering amplitude squared can be 
approximated as 



(2.7) 



spins 



Changing variables from t to u in (2.2), and using the 
bounds u m in = — s ~ — s and u max = — M 2 s/s ~ — M 2 , 
(2.7) contributes to dE/dx as 
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where we used 



I" n B (k) 
Jk 2fc 
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(2.8) 



(2.9) 



Adding (2.5) and (2.8) we obtain the muon energy loss 
to logarithmic accuracy, 



dEM e 4 T 2 r 



dx 



48tt 



ET 1 ET 



(2.10) 



2. Heavy quark energy loss in a QGP 

The heavy quark collisional energy loss arises from 
elastic scattering off thermal quarks (Fig. la) and glu- 
ons (Figs. lb-d). Compared to the QED case, there is 
one additional amplitude (Fig. Id), corresponding to t- 
channel exchange off thermal gluons. Using the cross sec- 
tions [20] corresponding to the QCD amplitudes of Fig. 1 
(with tree level gluon propagators in the ^-channel ampli- 
tudes), we can easily identify the origin of a logarithmic 
enhancement in the QCD analog of (2.2). There is a soft 
logarithm ~ J dt/t arising from |^Vf q | 2 and |.M g | 2 , and 
a collinear logarithm <~ J du/u from |Alg| 2 . The inter- 
ference terms are not enhanced by any logarithm - they 
however contribute to the constant to be evaluated in 
section II B. 



The heavy quark energy loss is readily derived from 
the QED result as follows. The contribution from |A4 q | 2 
is obtained by multiplying (2.5) by the number of quark 
flavours n/ and by the color factor {N 2 — l)/(4iV c ) = 
2/3. To get the ^-channel contribution from |.M g | 2 for 
scattering off thermal gluons, we also start from (2.5) and 
multiply by 1/2 since contrary to the electron, a gluon 
is its own antiparticle. This factor is compensated by a 
factor 2 arising from the difference between bosons and 
fermions when performing the integral over k (compare 
(2.6) and (2.9)). The color factor for this contribution 
is {N 2 — l)/2 = 4. Finally, the w-channel contribution 
from |A^g| 2 is obtained from (2.8) by multiplying by the 
color factor C 2 F = 16/9. Introducing the QCD coupling 



by e 



47ra s , we obtain the heavy quark energy 



loss in a QGP by summing all contributions, 



dE 
dx 



iiraiT 2 



/ n f \ET 2 ET 
1 + — In + - In — - + c(n f ) 
V 6 / m 2 n 9 M 2 y 11 



(2.11) 

47ra s T 2 (l + n//6) [18] is the Debye mass 
squared in the QGP. The constant c(n/) <~ O (1) is eval- 
uated in the next section, see (2.17). 



Here m 2 D 



B. The constant beyond leading logarithms 

In the QED case, the constant ~ O (1) in (2.10) was 
determined in Rcf. [17]. Here we infer from the QED 
calculation the constant c(nf) appearing in the QCD ex- 
pression (2.11). 

The logarithms \n(ET/m 2 D ) and ln(£T/M 2 ) arise 
from the ranges m 2 D -c \t\ <C s and M 2 <C \u\ <C s, 
with s ~ ET. This is why they could be easily ob- 
tained in section II A, using approximate expressions for 
the squared amplitudes (see (2.4) and (2.7)) in those 
kinematical domains. Controlling the constant next to 
the leading logarithms, however, requires considering 
the complete phase space < \t\ < |i m i n | — s and 
|u max | — M 2 < \u\ < |u min | ~ s. 

In order to treat correctly the domain \t\ ~ m 2 D , it is 
convenient to introduce an intermediate scale t* chosen 
as m 2 D < \t*\ < T 2 [17]. The contribution to dE/dx 
from \t\ < \t*\ is determined by the HTL self-energy of 
the exchanged gluon in Figs, la and Id, similarly to the 
QED case where it depends on the photon HTL self- 
energy (and where only Fig. la contributes). In both 
cases, the HTL self-energies have the same form, up to 
the replacement of the QED Debye mass by its QCD 
counterpart. Introducing the overall color factor Cf = 
4/3, the QCD result is thus inferred from [17] to be 



dE 
dx 



a s m 2 D 



\*\<\**\ 



In 



1**1 



In 2 



m 



D 



(2.12) 



We stress that the latter equation is valid beyond loga- 
rithmic accuracy, i.e. the constant next to the leading 
logarithm (written as In 2 here) is meaningful. 
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The contribution from \t\ > \t*\ can be evaluated by 
using tree level internal propagators in the amplitudes of 
Fig. 1. (Since \u\ > M 2 > T 2 , HTL corrections to the 
internal quark propagator in Figs, lb and lc are irrele- 
vant.) However, in order to control the constant beyond 
the leading logarithms, accurate expressions for the tree 
level cross sections have to be used. 

For the contribution from scattering off light quarks 
(Fig. la), the calculation is similar to the QED case, 
which is done by using 



oc 



t 2 



(2.13) 



instead of (2.4), and by performing the integrals in (2.2) 
with < \t\ < s. The QCD result reads: 



dE q 
dx 



4ira 2 T 2 n f 



\t\>\t*\ 



In 



8ET 



+ c 



, (2.14) 



where c = C'(2)/<(2) -7 ~ -1.147, with 7 ~ 0.577 being 
Euler's constant. 

Similarly, scattering off gluons (Figs, lb-d) should in 
principle be evaluated with the exact tree level cross sec- 
tion [20]. However, in the limit s ~ ET 3> M 2 we are 
considering, this cross section can be approximated as 



oc 



"s 2 
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1" 
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— u 
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J 2 




h 2_ 


+ 9 
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—u 



(2.15) 



The contribution of the first term of (2.15) to dE/dx is 
evaluated as for scattering off quarks (see (2.13)), except 
for the factor ns(fc) instead of tif (k) in (2.2). Up to the 
color factor 2/9, the second term of (2.15) has the same 
form as QED Compton scattering, and its contribution 
to dE/dx can be directly obtained from [17]. Summing 
the two contributions we get 



dEg 
dx 



AnatT 2 



\t\>\t*\ 



In 



AET 



2 / AET 
+ 9l ln ^ 



(2.16) 



Adding now the contributions from (2.12), (2.14) and 
(2.16) we obtain for the constant c(n/) in (2.11) 



c(n/) = an f + b ~ 0.146 n f + 0.050, 



(2.17) 



where the exact values of a and b are a = (2/3) In 2 
1/8 + c/6 and b = (31/9) In 2 - 101/108 + llc/9. 



III. IMPLEMENTING RUNNING COUPLING 

Implementing the running of a s in the calculation of 
dE/dx has already been done in Ref. [21] for the logarith- 
mically enhanced contribution from i-channel exchange. 
Here we will generalize the argument to the collinear log- 
arithm stemming from the it-channel. We also discuss 



the influence of running on the term oc a 2 c(n/) in (2.11), 
which allows us to define precisely the level of accuracy 
of our final result, see (3.7). 

We first note that the QCD result (2.11) (together with 
(2.17)) obtained with a fixed coupling suffers from a lack 
of predictability. Indeed, the fixed coupling calculation 
does not specify at which scale a s should be evaluated. 
As in PQCD calculations at zero temperature, 'fixing the 
scale' in a s demands to calculate the next order in the 
perturbative series for the observable of interest. This 
next-to-leading order (NLO) generally brings large log- 
arithms, which can be 'absorbed', via renormalization, 
by setting the scale of a s in the leading order result. In 
the case of i-channel scattering (Figs, la and Id), the 
vacuum contributions to the self-energy and vertex cor- 
rections conspire to yield a logarithmic dependence on 
the invariant transfer t. In order to avoid a large NLO 
contribution, one must set the scale of a s to ~ 0(|i|) 
[24]. 

In order to obtain a predictive result for the i-channcl 
contribution to dE/dx, we consider the term derived in 
section II A, which is logarithmically enhanced with a 
fixed coupling, and recalculate it with a running coupling 
a s (|i|). The procedure is rather trivial, 



ET 



<m 
\t\ 



£T ^(Ki). 
,2 t 



Using 



a. 



; (K|)=[4^ ln(|t|/A 2 )]- 



(3.1) 



(3.2) 



where O = (11 - jn f )/(An) 2 is the leading coefficient of 
the QCD /3-function, we see that using a running coupling 
amounts to perform the replacement 



a. In ■ 



ET 



m. 



a s (m 2 D )a s (ET) In 



ET 



m 



(3.3) 



D 



As noted in [21], the latter result becomes ^-independent 
in the E — > 00 limit, where the logarithmic enhancement 
of the fixed coupling result (2.11) is invalid. 

The accuracy of the above procedure is inferred by 
noting that the scale of a s in the r.h.s. of (3.1) can be 
in principle chosen as C|i| rather than \t\, with C a con- 
stant of order unity. This implies an ambiguity of order 
al In ET/m 2 D to the r.h.s. of (3.3) [25]. Thus specifying 
the term oc a 2 c(nf) in (2.11) is meaningful provided we 
assume 



ET 



ET 



a s (m 2 D ) In < 1 ln^-<ln^. (3.4) 



in 



D 



For the u-channel contribution (Fig. lc) to the cross 
section, a s should be evaluated at a scale ~ O (\u\). Sim- 
ilarly to the above discussion, the logarithmic integral 
appearing in (2.8) is modified to <~ J d|u|a 2 (|?i|)/|{(|, and 
the collinear logarithm in (2.11) becomes 



a„ In 



ET 

M 2 



a s (M 2 )a s (ET) In 



ET 
M 2 



(3.5) 
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Since M 3> mo and we already assumed (3.4), the uncer- 
tainty when choosing \u\ as the scale in a s is of relative 
order a s (M 2 ) In ET/M 2 <C 1 compared to the contribu- 
tion oc a 2 c(nf) in (2.11). 

After we have shown how the leading logarithmic terms 
in the fixed coupling result (2.11) arc modified when 
the running is taken into account (see (3.3) and (3.5)), 
we now discuss the modification of the constant term 



oc a 2 c(rif). Since this term is the difference between the 
full result (2.11) and the leading logarithmic terms, it 
should be clear from section II B that it is determined 
by integrals over t dominated by either \t\ <~ m 2 D or 
\t\ <~ ET, and by integrals over u dominated by |u| <~ ET 
[26]. There is no logarithmic spread there. Thus with 
running coupling, the term a 2 c(nj) should be evaluated 
at a scale chosen arbitrarily [27] between m 2 D and ET. 



We can now give the result for dE/dx of a fast heavy quark, at leading order in the running coupling, 



dE _ 4vrT 2 
dx 3 



ET 

1+ 'J±\ a s (m 2 D )a s (ET) In — + ^a s (M 2 )a s (ET) ln^- + a 2 c(n f ) 



ET 



m 



D 



M 2 



(3.6) 



This result has been derived for E 3> M 2 /T and assuming (3.4). The constant c(n/) is given in (2.17). As discussed 
above, we are free to set the scale in a 2 c(nj) as a 2 — > a s (m 2 D )a s {ET), in order to rewrite (3.6) as 



dE AttT 2 



-a s (m 2 D )a s (ET) 



dx 3 

where we displayed the order of neglected terms 



ET 2a s (M 2 ) l ET , N / , , N1 ET\ 
7 In— + c(n f )+0 (a s {m 2 D )\n—\ 

V m D / 



\ ^ 6 J m 2 D 9 a 8 {m 2 D ) M 2 



r 



(3.7) 



As a final remark, we note that when the logarithmic 
term oc In ET /m 2 D from i-channel exchange is dominant 
in (3.7), we obtain an interesting relation between the fast 
heavy quark collisional loss and the gluon Dcbyc mass, 



dE m 2 D ET 
-^—a s (ET)ln-^. 



(3.8) 



We used the self-consistent equation for the QCD Debye 
mass derived with running coupling [23] , 

m 2 D =47r(l + ^) as (m 2 5 )T 2 . (3.9) 

Using (3.2), the E — > oo limit of (3.8) yields an upper 
bound for dE/dx, 



dE 



< 



m, 



47T 



dx ~ 12tt/3 33 - 2n f 



m 



D ■ 



(3.10) 



For n f = 3 and m D ~ 0.66 GeV (at T = 0.2 GeV, see 
below), this gives the bound dE/dx < l.OGeV/fm. 



IV. A NUMERICAL ESTIMATE 



a s (m 2 D ) <C 1 in the (very) high temperature limit. The 
explicit terms in (3.7) are then all meaningful, including 
that oc c(rif). 

In practice, say under RHIC conditions, the values of 
m 2 D and ET are not very different (on a logarithmic scale) 
from A 2 . Therefore, the condition (3.4) is not satisfied in 
the strict sense, and the theoretical uncertainty affecting 

(3.7) is rather large. We take n/ = 3 and the generic val- 
ues (in GeV) A = 0.2, T = 0.2, M = 1.3 (charm quark), 
E = 20. We use (3.2) and (3.9) and find m D ~ 0.66 GeV 
(in good agreement with lattice QCD calculations, cf. 
[23]), a s (m 2 D ) ~ 0.58 and dE/dx ~ 0.6GeV/fm. The 
separate contributions in the bracket of Eq. (3.7) read 
[3.31 + 0.12 + 0.49], to be compared to the magnitude of 
neglected terms a s (m 2 D ) In ET/m 2 D ~ 1.29. 

Thus the theoretical uncertainty is larger than the 
collinear logarithm from u-channel exchange and than 
the constant c(n/) in (3.7) - making the approximation 

(3.8) reasonable under RHIC conditions. In the absence 
of an explicit calculation of the neglected terms, our re- 
sult dE/dx ~ 0.6GeV/fm in the above conditions might 
be accurate only up to a factor ~ 2 or so, and should be 
considered at best as a sound estimate. 



Our final result (3.7) for the fast heavy quark colli- 
sional energy loss is predictive: the scales at which to 
evaluate the different factors of a s are determined, and 
the order of its theoretical uncertainty is known. Con- 
sider the condition (3.4), in the limit of very high en- 
ergy E and very high temperature T, such that the log- 
arithm \aET/m 2 D <~ ln_E/T is kept fixed. Then the ne- 
glected contributions in (3.7) are indeed small, due to 
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